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In this paper we present the measurement of the integrated and differential production cross 
sections of pairs of prompt isolated photons in proton-proton collisions at a centre-of-mass energy 
of 7 TeV with the Compact Muon SoUenoid (CMS) detector at the Large Hadron Collider (LHC). 
A next-to-leading order (NLO) perturbative QCD (pQCD) prediction was compared to the mea- 
surement, which was performed on a sample corresponding to an integrated luminosity of 36 pb~^. 
While the agreement on integrated cross sections is satisfactory, a discrepancy is observed in the 
region of the phase space populated by photons with small relative angle, where the theoretical 
^ ^ prediction underestimates the measured cross section. 

in 



I. INTRODUCTION 



Studying the production cross section of pairs of energetic isolated photons at the LHC allows to test and 
^ validate pQCD calculations up to scales never probed before. Pairs of prompt photons, produced in the hard 
scattering of quarks and gluons, as opposed to photons produced in the decay of hadrons, also constitute the 
JL largest irreducible background to searches for a low mass Higgs boson decaying to two photons. 

D . 

In this paper we present the measurement of the integrated and differential production cross sections of 
'— ' pairs of prompt isolated photons performed with the CMS detector at the LHC, at a centre-of-mass energy 
^ of 7 TeV using the data collected in 2010 corresponding to an integrated luminosity of 36 pb~^. The cross 
^ sections are measured for two isolated photons with transverse energy greater than 23 and 20 GeV respectively, 
separated by more than R = 0.45 in the (77, 93)-plane and entering the pseudorapidity acceptance \ri\ < 1.44 or 
T— I 1-57 < I77I < 2.5. The isolation criterion is defined at parton level, we require less than 5 GeV of transverse 
energy to be deposited in a cone of radius R = 0.4 in the (77, </3)-plane. The differential cross sections are 
measured versus four variables, chosen either because of their sensitivity to the underlying QCD processes 
or their ability to describe the irreducible background in Higgs boson searches: the invariant mass m^^, the 
transverse momentum of the photon pair pT,'Y'y, the azimuthal angle difference between the two photons Atp^j 
and the cosine of the scattering angle in the Collins-Soper frame cos 9* [2]. 



^ After describing the event selection in section [H] we describe in section HI a novel method to discriminate 

*k>l prompt an non-prompt photons and give the details of the cross-section measurement in section [rv] The system- 
atic uncertainties are summarised in section [V] and the results of a NLO prediction compared to the measured 
^ cross section in section [Vl} Pairs of photons are called diphotons. 



II. EVENT SELECTION 



The CMS electromagnetic calorimeter (ECAL) consists of a cylindrical barrel part closed by two disc-shaped 
parts on either side, called endcaps. It is made of about 76000 lead tungstate crystals, that are slightly 
off-pointing with respect to the interaction point. Photons produced in CMS deposit their energy in ECAL, 
photon candidates are reconstructed by building clusters of ECAL crystals [3] . 

Events with two photons having transverse energies greater than 23 and 20 GeV respectively, separated by 
more than R = 0.45 and the clusters of which have \r]\ < 1.44 (ECAL barrel) or 1.57 < I77I < 2.5 (ECAL endcap) 
are selected. In order to reject non-prompt photons produced in jets, from the decay of energetic neutral mesons 
for instance, both photons are required to satisfy isolation criteria on the total track momentum and the total 
energy measured by the hadronic calorimeter (HCAL) in a cone around the photon candidate as well as criteria 
on the cluster shape. Furthermore it is required that no impinging track, defined as a track from the hard 
interaction with transverse momentum greater than 3 GeV, hits ECAL at less than R = 0.4 of the photon. In 
order to reject isolated electrons misidentified as photons, events with a hit in the first layer of the tracking 
detector matched to the photon cluster are rejected; the remaining contamination of the sample with Drell-Yan 
events is estimated from simulated events and corrected. 
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An acceptance correction factor is computed in each bin of the measured spectra, to account for the the finite 
detector resolution on the transverse energy and pseudorapidity, and the photon reconstruction and identification 
efficiencies. The efhciency of the isolation and cluster-shape criteria have been determined from data, using the 
tag and probe method [4 on electrons from Z boson decays. The efficiency of the selection on impinging tracks 
was determined with random cones, defined as a cone oi R = 0.4 around a direction at the same rj as the photon 
candidate and at a random in a 7r/2 window around the axis perpendicular to the photon direction, since the 
probability for an impinging track to hit a random cone is the same as the probability for an impinging track 
to hit the isolation surface of prompt photon. The acceptance correction factor is 76.2 ± 3.3% in the region 
corresponding to the entire pseudorapidity range. 



III. SIGNAL AND BACKGROUND DISCRIMINATION 



A sample of diphoton candidates selected with the requirements described in section |n] still contains a 
significant amount of background, that is diphoton candidates with one or both photons being non-prompt. In 
order to statistically determine the number of prompt isolated diphotons, we use the discrimination power of 
ECAL isolation variable X, defined as the sum of the transverse energies of crystals with Et > 300 MeV in an 
area around the photon ECAL impact position of outer radius 0.4, inner radius 3.5 crystal widths. Crystals 
belonging to the photon cluster or falling within a 5 crystal strip along Lp are removed from the sum. 

Given that prompt isolated photons are not produced along electromagnetic particles, and given that the strip 
along Lp excludes the deposits of potential conversion products, they have much smaller ECAL isolation values 
than non-prompt photons, that are produced in jets along other electromagnetic particles and the deposits 
of which tend to leak into the isolation region by multiple conversion and bremsstrahlung processes. The 
distributions for prompt an non-prompt photons of the variable I are significantly different and can therefore 
be used in a maximum likelihood fit to extract the number of prompt diphotons event in any given sam- 
ple of diphoton candidates, as shown in sectionjlVj The extraction of the two distributions is entirely data-driven. 
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FIG. 1: Validation of the T shape extraction on data for prompt (left) and non-prompt (right) photons. 
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The only contributions to the ECAL isolation of prompt photons come from pilcup activity and the 
underlying events. These contributions being independent of the distribution of the variable T of prompt 
photons can be extracted with the random cones defined in section |Tlj This random cone technique is validated 
with simulated diphoton events, as well as with data electrons from Z and W-boson decays, having radiated a 
small fraction of their initial energy in bremsstrahlung from interactions with the tracker and therefore leaving 
deposits in ECAL similar to those of prompt isolated photons. For W ev events, the sPlot technique was 
exploited, using a variable based on the missing transverse energy as discriminator. The comparison of the I 
distributions of W-decay electrons and random cones is shown on figure [T] 

The distribution of the variable I for non-prompt photons is extracted using a control sample containing less 
than 0.1 % of prompt photons, obtained from requiring one and one only impinging track to hit the isolation 
surface of the photon candidate. Subsequently the energy deposited by the charged particle that left the 
impinging track in a cone of radius 0.05 around its ECAL impact point is removed from the variable I which is 
then scaled to the nominal isolation area. This impinging track method is validated with simulated non-prompt 
photon events from QCD processes and with photon candidates in data, selected to be accompanied by two 
impinging tracks. The comparison of the I distributions of photon candidates with one impinging track and pho- 
ton candidates with two impinging tracks, the deposits of one of which have been removed, are shown on figure[l] 

Thus, using the data-driven random cone and impinging track methods, it is possible to extract the distri- 
butions of the variable I for prompt and non-prompt photons, in the ECAL barrel and the endcaps separately. 
The discrepancies observed in the validation procedure (Fig. [I]) will be taken into account in the systematic 
uncertainties. 



IV. CROSS-SECTION MEASUREMENT 



In each bin of the measured spectra, the number of prompt isolated diphotons is extracted by means of an 
extended maximum likelihood fit where the two-dimensional probability density functions (pdf) are defined 



as the product of the one-dimensional distributions /(X) extracted with the methods described in section III 



for each of the three types of diphotons: two prompt photons (the signal), two non-prompt photons and one 
prompt and one non-prompt photon (the background), . The likelihood function which is maximised is given 
in equation ([!]) where T are the three event types. A'' is the the size of the sample, Nt the numbers of events 
estimated in the fit for each type t, Ntot the sum over the three event types and ft{T^,I'^) is the probability 
for the ECAL isolation variables of the two photons to have values and for the given event type t. 



i=i teT 

Since the distributions of the variable I are extracted in the ECAL barrel and the endcap separately, the 
fit is performed for three different categories: for events with both photons in the ECAL barrel, events with 
both photons in the ECAL endcap and events with one photon in each part of the ECAL. As the distributions 
/(X) show a moderate dependence on rj and the pile- up and the distribution /(X) of non-prompt photons also 
depends on the transverse energy Et of the candidate, the events in the samples used for the extractions of 
/(I) are weighted so as to reproduce the distributions of 77, Et and the number of pile-up vertices of the fitted 
diphoton sample. The fit performed in the bin corresponding to 100 GeV < m^y < 140 GeV for both photons 
in the barrel is shown on figure [2] 

The number of prompt isolated diphotons extracted with the maximum likelihood method is known to be a 
biased estimator of the actual number of signal events, for fitted samples with small numbers of events. This 
bias is estimated with Monte-Carlo pseudo-experiments, found to be smaller than half the statistical error in all 
bins and negligible in the bins with more than 100 events, and the result of the fit corrected for it. The number 
of signal events estimated with the fit has been corrected for the reconstruction and identification efficiencies via 
the acceptance correction factor described in section |TT] and for the detector resolution by inverting the response 
matrix obtained from simulated events for m^^, pr.-yy, ^^-yy £^nd cos 9*. Given the excellent performance of 
ECAL, this matrix is nearly diagonal and no regularisation has been applied in the unfolding procedure. 
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FIG. 2: Projections of the two-dimensional fit in the bin corresponding to 100 GeV < m-^-^ < 140 GeV for photons with 
\ri\ < 1.44. The continuous blue hue represents the sum of the signal and the background contributions the dashed red 
line represents the background contribution only. In this bin, 161 diphoton candidates were selected and the number of 
signal events has been determined to be 72 ± 14. 



V. SYSTEMATIC UNCERTAINTIES 



The imperfect knowledge of the signal and background pdfs used in the fit is the main source of systematic 
uncertainties, which are estimated using Monte Carlo pseudo-experiments with varied distributions. The 
amplitudes of the variation is taken to be the discrepancies between the distributions extracted with the 
methods described in section III and the cross-check distributions, examples of which are shown on figure [l] 



They are of the order of ±0.01 for the signal and range from ±0.03 to ±0.05 for the background. 

The ECAL energy scale, known to 0.6 % in the barrel and 1.5 % in the endcaps [7J, affecting the definition of 
the acceptance via the Et thresholds defined in section |ll] and induced bin-to-bin migrations in the differential 
spectra affected by the energy measurement, that is m^^ and Pt, ■y-y- The former uncertainty is particularly 
important in regions populated by photons close to the Et thresholds of 20 and 23 GeV. 

The uncertainties associated with the computation of the acceptance correction factor, defined in section [ill 
are a combination of the statistical uncertainty on the samples used and the systematic uncertainties of the 
methods employed, taken to be the differences obtained when applying the method to samples of data and 
simulated events. The uncertainty on the integrated cross section is 3.3 % for diphotons entering the entire 
pseudorapidity acceptance. 

On top of these uncertainties, an uncertainty of 4 % is assigned to the knowledge of the integrated luminosity 
of the dataset [S] . The systematic uncertainties on the four spectra are summarised on figure [s] showing 
that with the exception of the region corresponding to m^j < 40 CeV, the measurement is dominated by 
statistical errors. The systematic uncertainties are dominated by the uncertainty on the signal and background 
pdfs, except for the same mass region, populated by photons close to the Et threshold, leading to significant 
uncertainties on the acceptance definition via the knowledge of the ECAL energy scale. 
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FIG. 3: Summary of the statistical and systematic uncertainties in the four measured spectra. 



VI. RESULTS 



A fixed order theoretical prediction is compared with the measured spectra. It includes all contributing 
processes at NLO: the quark annihilation process and the single and double fragmentation processes at order 
asO^ calculated with the diphox parton generator code ^ and the gluon fusion process, enhanced at the LHC 
because of high gluon densities, up to and including the order a^a^ computed with the GAMMA2mc generator 
code [TO]. 

As particles resulting from underlying event activity and hadronisation are not included in partonic event 
generators such as diphox and GAMMA2mc, the fraction of diphotons not selected due to underlying hadronic 
activity falling inside the isolation cone has been estimated using the pythia 6.4 [TT] event generator. This 
results in the parton-level cross section being reduced by 4.7 ± 0.3 %. Theoretical uncertainties on the 
knowledge of the parton distribution functions (PDF) and as were determined following the PDF4LHC 
recommendations |12| and the factorisation, renormalisation and fragmentation scales were varied by factors 
1/2 and 2, keeping the ratio between two scales less than 2. The predicted integrated cross section is: 
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a{pp 77)1 



prcd. 

|»)|<2.50 ~ 



52.7 



+5.8 
-4.2 



(scales) ±2.0 (PDF) pb. 



(2) 



It is found to be compatible with the measured integrated cross section within the statistical and systematic 
uncertainties: 



cr(pp 77)1 



mcas . 
|r;|<2.50 



62.4 ± 3.6 (stat) 1^% (syst) ± 2.5 (lumi) pb. 



(3) 



The predicted differential cross-section spectra for the variables my^, pt,77, ^^-yy and cos 0* are compared 
to the measured spectra on figures |4j [5j [6] and [7| respectively. 



Figure [6] shows that the theoretical cross section is underestimated in the region Aip. 



77 



< 2. 



Indeed, in 



the leading order (LO) description of the gluon fusion and quark annihilation processes, the two photons are 
back-to-back because of momentum conservation. Therefore the LO term does not contribute to this phase 
space region, which is effectively covered in the NLO calculation by one order only. This affects in particular 
the fragmentation contributions [13 . Diphotons with small relative angle between the photons populate the 
regions m^^ < 40 GeV and pt,77 > 40 GeV (because of the requirement of Et > 20, 23 GeV) both of which 
are also underestimated by the theoretical prediction. 

Analogous disagreements were observed in the diphoton production in proton-antiproton collisions by the 
CDF [H] and DO jl5J collaborations as well as in proton-proton collisions by the ATLAS 16 collaboration. 
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FIG. 4: Measured cross section of diphoton production as a function of the invariant mass of the photon pair (left) and 
bin-by-bin comparison with the theory (right) for photons within the pseudorapidity region \ri\ < 1.44 or 1.57 < |r;| < 2.5. 
The total systematic uncertainties are represented by the shaded area, the different contributions are added in quadrature 
sequentially. 



VII. CONCLUSION 



The integrated and differential production cross sections of isolated photon pairs have been measured in 
proton-proton collisions at a centre-of-mass energy of 7 TeV. The differential cross sections were measured as 
functions of the diphoton invariant mass and transverse momentum, the azimuthal angle difference between 
the two photons, and the cosine of the scattering angle in the Collins-Soper reference frame. The non-prompt 
photon background contamination from hadron decay products has been estimated with a statistical method 
based on an electromagnetic energy isolation variable, the signal and background distributions of which have 
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FIG. 5: Measured cross section of diphoton production as a function of the transverse momentum of the photon pair (left) 
and bin-by-bin comparison with the theory (right) for photons within the pseudorapidity region \rj\ < 1.44 or 1.57 < 
|?7| < 2.5. The total systematic uncertainties are represented by the shaded area, the different contributions are added 
in quadrature sequentially. 
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FIG. 6: Measured cross section of diphoton production as a function of the azimuthal angle between the two photons (left) 
and bin-by-bin comparison with the theory (right) for photons within the pseudorapidity region I77I < 1.44 or 1.57 < 
|?7| < 2.5. The total systematic uncertainties are represented by the shaded area, the different contributions are added 
in quadrature sequentially. 



been entirely extracted from data resulting in systematic uncertainties of approximately 10 % on the measured 
cross section. The measurements have been compared to a theoretical prediction performed at NLO order 
accuracy using the state-of-the-art fixed order computations, showing that in regions of phase space mostly 
populated by photons emitted at small relative angle the prediction underestimates the measured cross section. 
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FIG. 7: Measured cross section of diphotori production as a function of the absolute value of the cosine of the polar angle 
between the two photons in the CoUins-Soper reference frame and the photon pair (left) and bin-by-bin comparison with 
the theory (right) for photons within the pseudorapidity region I77I < 1.44 or 1.57 < I77I < 2.5. The total systematic 
uncertainties are represented by the shaded area, the different contributions are added in quadrature sequentially. 
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